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Abstract — In this paper, we investigate the problem of de- 
signing compact support interpolation kernels for a given class 
of signals. By using calculus of variations, we simplify the 
optimization problem from an infinite nonlinear problem to 
a finite dimensional linear case, and then find the optimum 
compact support function that best approximates a given filter 
in the least square sense (£% norm). The benefit of compact 
support interpolants is the low computational complexity in 
the interpolation process while the optimum compact support 
interpolant gaurantees the highest achivable Signal to Noise Ratio 
(SNR). Our simulation results confirm the superior performance 
of the proposed splines compared to other conventional compact 
support interpolants such as cubic spline. 

Index Terms — Spline, Interpolation, Filter Design 



I. Introduction 

DUE the existence of powerful digital tools, nowadays it is 
very common to convert the continuous time signals into 
the discrete form, and after processing the discrete signal, we 
can convert the discrete signal back to the original domain. 
The conversion of the continuous signal into the discerete 
domain is usually called the sampling process; the common 
form of sampling consists of taking samples directly from 
the cotinuous signal at equidistant time instants (uniform 
sampling). Although the samples are uniquely determined 
by the continuous function, there are infinite number of 
continuous signals which produce the same set of samples. 
The reconstrcution process is defined as selecting one of the 
infinite possibilities which satisfies certain constraints. For a 
given set of constraints, a proper sampling scheme is the one 
that establishes a one-to-one mapping between the discrete 
signals and the set of continuous fucntions that satisfy the 
constraints. One of the well-known constraints is the (finite 
support in fourier domain) condition [?]. Due to the advances 
in the wavelet theory [?], [?], [?], a new intrest in continuous- 
time modeling and spline interpolation has been generated. 
Multiresolution analysis [?], [?], self-similarity [?], [?], and 
singularity analysis [?] are inseparable from a continuous- 
time interpolation. Although proper sampling schemes provide 
a one-to-one mapping between the continuous and discerete 
forms, we still need the tools to recover a continuous time 
signal from its samples. In fact, interpolation techniques using 
splines are one of the best options here. 

In this field, polynomial splines, such as B-Splines, are 
particularly popular; mainly due to their simplicity, compact 
support, and excellent approximation capabilities compared to 
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Fig. 1. Sampling process modeled by multiplying an impulse train by a 
continuous time signal. 



other methods. B-Spline interpolations have spread to various 
applications [?], [?], [?]. 

Many advantages of the B-splines arise from the fact that 
they are compact support functions. However, there is no 
evidence that they are the best compact support kernels for the 
interpolation process; i.e., it may be possible to improve the 
performance without compromising the desired property of the 
compact supportedness. In this paper, we focus on the problem 
of designing compact support splines that best resemple a 
given filter such as the ideal lowpass filter; more precisely, we 
aim to find a compact support spline that minimizes the least 
squared error when its cardinal spline is compared to a given 
fucntion. The given filter may be any arbitrary function that 
reflects the properties and constraints of the class of signals 
that enter the sampling process. 

The remainder of the paper is organized as follows: The 
next section briefly describes the spline interpolation method. 
In section Hill a novel scheme is proposed to produce new 
optimized splines for interpolation regardless of the type of 
filtering. The performance of the proposed method is evaluated 
in section|lV]by comparing the interpolation results of the pro- 
posed method to those of well-known interpolation techniques. 
Section [V] concludes the paper. 

II. Preliminaries 

In this paper, the following notations and definitions are 
used: 

Definition 1. For a continuous-time signal x(t), a continuous- 
time signal x p (t) and a discrete-time signal Xd[n] are defined 
as follows: 

x d [n] = x(nT) (1) 

OO 

x p {t) 4 x{t)p{t) = Xd[n]5{t-nT) (2) 

n — — 00 

where p(t) = J2n=-oo — ^ s me periodic impulse train. 
The sampling period T is normalized to 1 without any loss of 
generality. 



The sampling process is shown in Fig [T] 
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Fig. 2. Spline of the degree 3 conditions. 

Definition 2. A Linear Time Invarient (LTI) filter with impulse 
response h(t) is called to have interpolation property if and 
only if, 

h p (t) = S(t). (3) 

Definition 3. For a continuous-time signal x(t) and any odd 
integer m, x™(t) is a polynomial spline of degree to if, 

1) for any n G Z, x™(t) is a polynomial of (at most) degree 
to in the interval [n, n + 1] (Compact support), 

2) for any n G Z, x" l (n) = Xd[n] (Interpolation property), 

3) x? G C ,m - 1 (-oo,oo) (Smoothness). 

According to the first property, (m + l) t/l derivative of x™ 
is equal to an impulse train. 

Definition 4. For the polynomial spline x™(t), the polynomial 
spline coefficients £™[n] are defined as: 

00 (j m + 1 
x™(t)= ]T ±?[n]6(t-n)±-~x?(t) (4) 



where (u p n+1 ) 1 (t) is defined as the inverse of u r p n+1 {t), i.e, 
* (t) = 5(t), and Xf(z) is the z-transform 

of x™[ra]. 

Proof: From Def. El 

= * i™) (i)p(i) 

= {u™ +1 *x™)(t) (7) 

Hence, 

The ROC of ?7™ +1 (z) is \z\ > 1 and there are no zeros in this 
region. Since the ROC of X'T (z) in not enclosed in the unit 
circle, (C/™ +1 )~ 1 (z) and X7 l (z) have a region in common. 
Thus, 

±™(t) = ((u™ +1 )- 1 *x™)(t) (9) 

and according to (O, 

a;™(t) = (u m+1 * x™) (t) 

= (10) 

■ 

Definition 5. A discrete-time signal yd[n] is called a proper 
signal if and only if it is bounded and has a unique and 
bounded inverse 1 [n] . 

Definition 6. For any continuous-time signal y(t), if yd[n] is 
a proper signal, then y(t) is defined as follows: 

y(t) = ((y P )~ 1 *v)(t) (id 

Corollary 1. y(t) in ( fTTT ) is the impulse response of a filter 
with interpolation property, in other words: 

y p (t) = 6(t) (12) 



To determine all polynomials of degree m that form x™(t), 
the m + 1 unknown spline coefficients should be found in 
order to satisfy the conditions 2 and 3 (Fig. |2}. If the goal is to 
discover a piecewise polynomial signal that is (to — 1) times 
differentiable with continuous derivatives, a natural method 
is to derive x™[n] according to Xd[n] and then calculate the 



integral of x™(t), m 
<(t) = 



1 times, i.e, 

■ti 



±™(t )dt ■ ..dtm-ldt r . 



* x r 



(t) 



(5) 



where u^it) is the unity step function and for any k G N, 

Lemma 1. If the Region Of Convergence (ROC) of X d n (z) is 
not enclosed in the unit circle (i.e, there exists z G ROC {X r d n } 
such that \z\ > 1), then from Def. [3] x™{£) will be uniquely 
determined from Xd[n], and 



'(*)=(( 



u m+1 * (u 



m+l\ 
P ) 



(t) 



Proof: 



y P (t) = y(t)p(t) 

= ((y P r 1 *y P )(t) = S(t) 



(13) 



Lemma 2. For a given proper signal y™[n], the signal y r s n {t) 
is only a function of m and does not depend on y r d n [n] . 

Proof: 



yT(t) 



(((yDp)- 1 *^) (*) 

((yj-^y?) (t) 

((j/p)- 1 *^ 1 *^ 1 )- 1 )*^) (t) 



( 



,m+l 



(14) 



(g) Hence y™(i) is only a function of to. 
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Definition 7. According to the above corollary, c m (t) = y™ (i) 
is defined as the cardinal spline of degree m (Fig. |3)- 

From © and (TBI it can be concluded that the polynomial 
spline interpolation is a linear shift invariant process according 
to Xd[n] and its impulse response is represented by c m (t). 
c m (t), on the other hand, can be derived according to any 
arbitrary polynomial spline y™(t) such that y™(t) is a proper 
signal, i.e, 

xT(t) = (c m *x p )(t) 

= (l/r*((l/p) _1 *a:p))(t) (15) 

The above equation divides the whole interpolation process 
into a discrete-time and a continuous-time section. If is 
chosen as a time limited basis, both sections of this process 
can be extremely simplified and a considersble amount of 
continuous-time calculation can be avoided. 

Theorem 1. k = m + 1 is the least positive integer for which 
there exists a polynomial spline of degree m, such as y™(t), 
that vanishes outside the interval (0, k), i.e, 



(16) 



then k = m + 1, 



Proof: Suppose that y" l (t) satisfies (O and Y™(z) is the 
z-transform of y™[n], where y™[n] is the coefficients signal 
for the polynomial spline y" l (t) according to the definition 
©. It can be claimed that, 



(17) 



In order to prove (TTTt . we define the sequence of polynomi- 
als {QJ™ o such that Q (z) = Y^iz' 1 ) and for 1 < n < m, 



Qi = z 



dz 



£ 

71 = 



y?[n]n l z n 



(18) 



Also polynomial H is defined as follows: 




m! i! — ' 

i=0 

m / ^ 

i__o \n= 
' n=0 i=0 ^ ' 

= i£*n»](i-»r d9) 



Thus, according to (0, for any t > k, H(t) is equal to y™(t), 
i.e, 

Vt> k^ H{t) ^ y™(t) = (20) 

Since iJ is a polynomial that vanishes for infinite values of £, 
it should be equivalent to zero: 



H(t) = 



,(l) = 0m- l (l) 



Q (1) = (21) 



From the above equations, it is concluded directly by induction 
that, 



Vn e N; < n < m 



dz' 



-YT^- 1 ) 







(22) 



Hence, (z — l) m+1 divides YJ™^ 1 ). On the other hand, since 

Y^i^ 1 ) = ELo shows that + 1] ^ o. 

Hence, 

G(m+l-e,m+l+e) ^ (23) 

Thus k > rn + 1. 

Finally, it must be shown that there exists a polynomial 
spline of degree m that is supported in the interval (0, m + 1). 
Suppose that F d m (2) = - l) m+1 , then 



V?[n] = (-l) n 



rn + 1 

n 
m+l 



y?(t) = u 



— „, m + i 



m+l 



C(*)= ^y7[n}u m+1 (t-n). 



(24) 
(25) 

(26) 



Thus, for all t > m + 1, =0. ■ 

Definition 8. The polynomial B-Spline of degree m is defined 
as follows: 

m+l , | i \ 

nt) ^E(-i)- r « m+1 (t - n) (27) 



71=0 



Now according to the lemma [2] the polynomial spline in- 
terpolation process can be implemented by a compact support 
kernel such as the polynomial B-Splines, i.e, 

(28) 



(29) 



oo 



n= — oo 



III. Proposed Optimized B-Spline 

In many applications, it is desirable that the interpolation 
filter resembles an ideal filter, and the second and third prop- 
erties in Def. [3] may be ignored. In this section an optimized 
basis spline will be introduced to emulate a desired filter. 

Definition 9. Let D denote the set of all continuous-time sig- 
nals that satisfy the dirichlet conditions, i.e, for any y(t) 6 D 

1) y{t) has a finite number of extrema in any given interval, 

2) y(t) has a finite number of discontinuities in any given 
interval, 

3) y(t) is absolutely integrable over a period, 

4) y(t) is bounded. 

First of all, an affine subspace of all signals that satisfy the 
dirichlet conditions will be defined, and then the optimized 
solution will be obtained in this set by the calculus of variation. 

Definition 10. Let yj{n\ be a proper signal that vanishes for all 
n < and n > m+l, then x m (yd) is the set of all continuous- 
time signals y(t) that satisfy the following conditions, 
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Fig. 3. Cardinal splines of different degrees. 



1) yeD 

2) Vn G N; y(n) = y d [n] 

3) Vi $ (0,m + l); y{t) =0 

The use of y(t) G X m (yd) as a basis spline to interpolate 
Xd[n] according to equations (l28l l and ( f29b is a linear time 
invariant process with the impulse response y(t). 



Definition 11. The error function e x : x m {Vd) 
as follows: 



is defined 



e x {y) 



y*x p -x\\ 2 

\(y*x p )(t)-x{t)\ 2 dt 

) 

/oc 
\F{y*x p }- T{x}\ 2 df 
-OO 

\T{{(y p )- 1 ^y)*x p }-T{x}\ 2 df 

) 

,T{x p } 



F{y P } 



F{y}-F{x}\ 2 df. 



(30) 



where T is defined as the continuous time Fourier transform 
operator. 

Definition 12. According to the above definition, if is a 
proper signal that vanishes for all n < and n > m + 1, an 
optimized basis spline p m [x,p™] is defined as follows: 



p m [x,p%] = argmin e x {y) 



(31) 



Now, for a given proper discrete signal with the required 
vanishing property, we employ the calculus of variations in 
order to find the optimized continuous basis spline p m that 
minimizes the error function e x (p m ). 

Theorem 2. Equation OTb has a unique solution that satisfies 
the following property, 

for all t e (0,m + 1), where y(t) = y*(-t). 

Proof: For 7 e X™(0) and any e > 0, we have p m +£7 G 
X m (/3™), and the variational derivation of e x (p m ) with respect 



to p" 1 with 7 as the test function is equal to 

(e a (p ),7) = lim 

OO s 

- 2 /™ r 

—00 v 

7 ! " 





The proof of (l33l is presented in (34). Since X n "{P'd) ls 
boundless, in order to minimize e x (p m ), (e(p m ),j) should 
be zero for all 7 G X m (0)> which implies that the second term 
inside the integrals should be zero for t G (0, rn + 1), i.e, 



^{Xp} 



Hp m } 
Hp?} 



T{x p } - F{x} 



(35) 



The above equation directly yields 



Thus, it is proven that the optimized basis spline which gives 
the best estimation of x, should satisfy (I32t . By defining 

(36) 
(37) 



w P (t) = (^*^)*[(P™)- I *(^ r )^ 1 ] 



1-1 



we can rewrite (l32l as (v p * p 



m+l) 

this equation is only valid in a particular interval, (v p ) 



)(t) = w(t)\ te m m+1 y Since 
-1 



cannot be used to obtain p m . However, since v p is an impulse 
train (convolution of four impulse trains) we will show that the 
continuous functional equation in (l32l l boils down to solving a 
finite Hermitian system of linear equations. For this purpose, 
we first define the two following sequences of functions which 
are supported only on a unit-length interval: 



Rn(t) = 
W n (t) = 

Now, here is (l32l in the matrix form 




(38) 
(39) 

(40) 



Since Vd{— n] = v* d [n], the above matrix V = [vd[i — 

j]]i=i....,m+i;j=i m+i is Hermitian also. Now by solving 

(l40l >. {R n }™ =0 is derived and thus the optimized basis spline 
is given as 



Mo] 


M-i] ■ 


v d [—m] 




r -Ro - 




r Wo -1 


Mi] 


MO] ■ 


■ v d [-m+l] 




Ri 






v d [m] 


M«i-i] • 


MO] 




-R m - 




-W m . 



p m {t) = Y J Rn(t-n) 



(41) 



ri=0 
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{e x (p m [x,Pd]),i) = lim 



Mp m +£7)-e4p m ) 



= lim - 

E^O £ 



= lim - 

e^o e 



/ P{p-+ £ 7} 
9? 



3? 



| ^{p m } 

V-Hp™} 

p{p m } 

\Hp?} 



H x p} — -^{x} d/ 



H x p} ~ H x } 



+ 



lim 



V-Hp?} 



4f 



lim - 

e->0 e 



+ 3 



= lim - 

e-s-0 £ 



+ 3 



- r I* f - 3*1) nxp>) * ( Ly +£7 > + y^ .Ffe} - 2^ { , } 



/ ^{ P - + e7 } T{p m } 

\H(P m + E7M -^{P? 1 } 
3? 

^^{p m + £ 7}--F{p m } 



-F{a;p} 



V -Hp™} 



H x p} \ 3? 



^FjfT -:- ;-,}-:- J ■{,,■■•} 



Hp? 1 } 

"{(«)*'!«{( 



■Hppn 
^{ P " i + e7 } + ^{p m } 



J^} - 2.F{a:} 



^{2p m + £7 } 



- 2T{x} }df 



Hp?} 



\Hp?} 



df 



I K { ^-^l * { 2 mF} Hxv} - 2Hx} ] + 3 {h^} Hxp} } *{ 2 mr} Hxv} - 2Hx} 



df 



= 23? 
= 23? 
= 23? 
= 2 



^{xp} 



•^{p m } 



7W-F- 1 
3?{7(*)}3?< J" 1 



-^"{Pp"} 

Hp m } 



T{x p } - T{x} 



Hp?} 



H x p} - H x } 



df 
df 



T{x P } 



Hp?} 
Hx P } 



H x p} 
Hp?} 

T{x p } 



3{ 7 «H T 



Hp?} 
H*p} 



Hp?} 



Hp m }-H*} 

T{ P m }-T{x} 

Hp m }-H*} 



dt 



Hp?} 
H x p} 



<it 



Hp?V 



dt 



(34) 



This optimized basis spline minimizes the mean squared 
error of interpolation. On the other hand, smoothness of the 
optimized basis spline and causality of the prefilter can be 
achieved by adjusting pd- 



Another application of (1321 is to approximate an ideal 
interpolation filter by an optimized basis spline. In fact, basis 
splines are superior to FIR filters. 

Now the goal is to design p m such that p rn would be the 
best estimation of h, which denotes the impulse response of 
a filter that has the interpolation property. 

Lemma 3. (Estimating a desired filter) Assume h(t) is the 
impulse response of a linear time invariant filter which satisfies 
the interpolation property and let p™[n] be a proper discrete 



signal, then, 



arg mm 



\h-y\\ 2 =P m [h,p7] 



Proof: 

eh(y) 



/OO 
H{y*h P }-Hh}\ 2 df 
-OO 

/oo 
\Hy}HK) - Hh}?df 
-oo 

/oo 
\Hy}-Hh}\ 2 df 
-oo 



-00 

\\h-y\ 



(42) 



(43) 
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Fig. 4. The optimized spline versus cubic B-spline. p s {h,p|} is the opti- 
mized basis spline built for estimating the ideal lowpass filter h(t) 



a(nt) 



with P^(z ) = Q-235Z + 0-484z 2 + 0-235.Z 3 




Fig. 5. The comparison of the proposed method with the cubic spline for 
an ideal lowpass filter. 



Thus 



argmin \\h - y\\ 2 = argmin e h (y) = p m [h 7 pf] (44) 



Corollary 2. For an impulse response hit) with the interpo- 
lation property, if p m G X m (p™) denotes the optimum basis 
spline for which p m best approximates h(t) (i.e., p m (t) = 
argmin \\h — y\\ 2 ), then p rn £ X m (Pc?) satisfies: 

The proof follows directly from (l3Zt and the fact that h(t) has 
the interpolation property. 

Figure |4] shows the cubic B-spline and the optimized B- 
spline designed for estimating the ideal lowpass filter h(t) = 
sinc(t) with p 3 {z) = 0.235z + 0.484z 2 + 0.235z 3 , while Fig. 
[5] shows p 3 [h, Pj]{t) in comparison to c 3 (t). 

IV. Simulation Results 

To compare the proposed method with the existing inter- 
polation techniques, we have performed various simulations. 
The cubic B-spline, due to its short time support and relatively 



high accuracy in approximating the ideal lowpass filter, is the 
most common technique for interpolating 1-D lowpass signals. 
For the purpose of comparison, we have optimized a spline 
function with the same time support as an ideal lowpass filter. 
Figures |4] and [5] show the shape of the obtained B-spline and 
the interpolating spline, respectively. Figure |4] shows that the 
energy is more concentrated in the middle of the cubic B- 
spline while the optimized spline has slower decaying rate 
of energy at the sides. The resultant interpolation splines, 
as depicted in Fig. [5] reveal that the main advantage of the 
optimized spline (p 3 ) compared to the cubic spline (c 3 ), is the 
smaller error in the first side-lobe. The SNR values of p 3 and 
c 3 with respect to the sine function are 20.39 and 13.15dBs, 
respectively. 

For a more realistic comparison, we have applied different 
interpolation techniques on standard test images. For this 
purpose, the original images, with or without applying the 
anti-aliasing filter (ideal lowpass filter), are down-sampled by a 
factor 2 in each direction (25% of the original pixels) and then 
they are enlarged (zooming) using the interpolation techniques. 
The comparison is made with the following interpolation 
methods: 1) bilinear interpolation, 2) bicubic interpolation, 3) 
wavelet-domain zero padding cycle-spinning [?], and 4) soft- 
decision estimation technique for adaptive image interpolation 
[?]. Also for our proposed method, two different scenarios 
are implemented: the ideal filter for which we are optimizing 
the spline function is first taken as a sine filter, and first 
the spectrum of the original image. In order to evaluate the 
quality of the interpolated images, we have considered the 
Peak Signal-to-Noise Ratio (PSNR) criterion. Table ^indicates 
the resultant PSNR values when the original image is subject 
to the anti-aliasing filter before down-sampling while the error 
is calculated based on the image without applying the filter. 
Table HI] contains similar values while the basis for the error 
calculation is the anti-aliased image. In both cases, the PSNR 
criterion favors the proposed optimized spline for the ideal 
lowpass (sine) filter. On the average, the proposed method out 
performs the other standard interpolating methods by 0.74dB 
in Table U and 4. 19dB in Table M 

To exclude the effect of the anti-aliasing filter, the simu- 
lations are repeated without applying it and the results are 
presented in TableHTTl As expected, the optimized spline which 
is matched to the spectrum of the original image outperforms 
other competitors in all cases. It should be mentioned that the 
function p m \x 7 p^\ is not a universal filter in this case and 
depends on the choice of the image. 

Although the PSNR value is a good measure of global 
quality of an image, it does not reflect the local properties. 
In order to present a qualitative view of various interpolation 
methods, we have plotted the enlarged images for a segment 
of the Lena test image in Fig. [6] To highlight the differences, 
one could compare the texture on the top and the sharpness 
on the bottom edge of the hat. 

V. Conclusion 

The interpolation problem using uniform knots is a well 
studied subject. In this paper, we considered the problem of 
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TABLE I 

PSNR (dB) Results of the Reconstructed Images by Various Methods, The Original Image was Anti-Aliased Before Sampling and 
The Results are Compared to The Original Image (Image Enlargement from 256 x 256 to 512 x 512) 



Images 


Bilinear 


Bicubic [?] 


WZP-CS [?] 


SAI [?] 


Optimized Spline for 
Image p 3 [x,p 3 ,] 


Optimized Spline for the Ideal 
LowPass p 3 [sinc(t) t p^] 


Lena 


30.33 


30.44 


30.12 


30.96 


32.46 


32.20 


Baboon 


22.40 


22.52 


22.41 


22.89 


22.20 


24.22 


Barbara 


24.39 


24.42 


24.34 


24.65 


24.33 


25.14 


Peppers 


29.46 


29.46 


29.14 


29.72 


31.14 


31.04 


Girl 


30.98 


30.98 


30.66 


30.43 


31.70 


30.87 


Fishing bout 


27.68 


27.48 


28.12 


30.92 


28.61 


29.71 


Couple 


27.35 


27.50 


27.27 


27.39 


28.04 


28.97 


Overall Average 


27.51 


27.54 


27.43 


28.14 


28.35 


28.88 



TABLE II 

PSNR (dB) Results of the Reconstructed Images by Various Methods, The Original Image was Anti-Aliased Before Sampling and 
The Results are Compared to The Anti-Aliased Image (Image Enlargement from 256 x 256 to 512 x 512) 



Images 


Bilinear 


Bicubic [?] 


WZP-CS [?] 


SAI [?] 


Optimized Spline for 
Image p 3 [x,p 3 d ] 


Optimized Spline for the Ideal 
LowPass p 3 [sinc(t), p 3 ] 


Lena 


31.56 


31.72 


31.62 


32.45 


35.28 


35.21 


Baboon 


26.88 


27.26 


26.89 


28.48 


26.33 


35.70 


Barbara 


30.61 


30.74 


30.65 


31.86 


30.40 


35.72 


Peppers 


31.62 


31.82 


31.77 


32.14 


37.41 


38.02 


Girl 


34.08 


34.24 


34.25 


33.10 


35.89 


34.28 


Fishing bout 


29.91 


30.19 


29.93 


30.92 


32.00 


34.87 


Couple 


29.81 


30.10 


29.84 


29.77 


31.44 


34.28 


Overall Average 


30.64 


30.87 


30.71 


31.25 


32.68 


35.44 



TABLE III 

PSNR (dB) Results of the Reconstructed Images by Various Methods, The Original Image was not Anti-Aliased and The Results 

Are Compared to The Original Image 256 x 256 to 512 x 512) 



Images 


Bilinear 


Bicubic [?] 


WZP-CS [?] 


SAI [?] 


Optimized Spline for 
Image p 3 [x, p 3 ] 


Optimized Spline for the Ideal 
LowPass p 3 [sinc(t) , p 3 ] 


Lena 


30.21 


30.13 


30.05 


30.88 


32.29 


30.95 


Baboon 


21.67 


21.34 


21.70 


22.09 


22.50 


21.63 


Barbara 


23.90 


23.32 


23.88 


23.71 


25.10 


22.58 


Peppers 


28.82 


28.61 


26.93 


28.91 


30.64 


29.77 


Girl 


30.41 


29.97 


30.20 


29.94 


30.90 


29.20 


Fishing bout 


27.10 


26.93 


27.07 


27.63 


28.50 


27.66 


Couple 


26.92 


26.73 


26.86 


26.93 


27.91 


27.08 


Overall Average 


27.00 


26.72 


26.67 


27.16 


29.12 


26.98 



optimizing the interpolation kernel for a given class of signals 
(represented by a filter). Although functional optimization 
in the continuous domain is often very difficult, we have 
demonstrated the equivalency of this problem with a finite 
dimensional linear problem which can be easily solved using 
linear algebra. As a special case, we considerd the class of 
lowpass signals which is associated with the sine function as 
the optimum interpolation kernel. For the optimum compact 
support interpolant, we compared our function with the con- 
ventional cubic B-Spline; the simulation results indicate ldB 
improvement in the SNR of the interpolated signal (on the 
average) using the introduced function, and 7dB improvement 
compared to the cardinal spline itself (compared to the sine 
function) (Fig. |5}. 
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(d) (e) (f) 

Fig. 6. Comparison of different methods for the Lena image: (a) The original image, (b) bilinear interpolation, (c) bicubic Interpolation, (d) WZP Cycle- 
Spinning [?], (e) SAI [?], and (f) the proposed method. 



